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Preface 

This  note  is  an  elaboration  of  the  third  note, 
"Transients,"  in  the  mimeographed  series  of  195i4-«  We 
consider  certain  matheinatically  simple  flows  which  are 
explicitly  solvable  and  are  chosen  to  exhibit  certain 
basic  properties  of  plasmas.   In  particular,  we  note 
that  gyro  and  plasma  oscillations  are,  in  the  simple 
Case  considered,  manifestations  of  a  single  phenomenon 
which  could  be  termed  "gyro-plasma"  oscillations. 
Elementary  properties  of  the  Debye  shielding  length, 
Larmor  radius,  and  another  quantity  termed  "magnetic 
shielding"  length  are  discussed. 

The  assistance  of  Albert  Blank,  Ignace  Kolodner, 
and  Paul  Reichel  in  preparing  this  note  is  gratefully 
acknowledged. 
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Special  Solutions 

In  this  note  we  consider  certain  special  flows  which 
are  simple  enough  to  be  obtained  explicitly  and  which 
illvuninate  certain  general  fluid  properties  connected  with 
the  terms  g3rro  (or  Larmor )  frequency,  plasma  frequency, 
and  Debye  length.   They  have  the  feature  in  common  of 
dependence  on  a  single  independent  variable,  either  time 
or  space.   In  particular  we  shall  conclude  that  gyro  and 
plasma  oscillations  are  different  facets  of  a  single 
phenomenon,  rather  than  distinct  phenomena. 
1.   Transients 

First,  let  us  consider  a  space-independent  problem 
with  time  alone  as  the  variable.   Prom  B  +  curl  E  =  0  we 
conclude  that  B  is  constant  in  time  also.   Prom  dlv  D  =  q 
we  conclude  that  the  medium  is  neutral.   There  remains 

(1)  KE  =  -  $:  q^  u^   ,   I-qj,  =  0 

as  the  only  non-trivial  one  of  Maxwell's  equations.   From 
the  conservation  of  mass,  we  conclude  that  p^  and  conse- 
quently  q^  are  constant  in  time.  W6  take  the  momentijin 
equations  with  frictlonal  terms  proportional  to  velocity, 

(2)  p^  Uj,  =  q^(E  +  u^  X  B)  -  j;  a^^  (u^  -  u^  ) 

s 

a   =  a„  >  0   . 
rs    sr 
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Since  p  ,  q  ,  and  B  are  constants,  (1)  and  (2)  are  a  set 
of  linear  ordinary  differential  equations  for  E  and  u 
as  functions  of  time. 

It  is  convenient  to  consider  first  the  frictionless 
case,  a  =  0«  Also,  as  a  formal  procedtu?e,  let  us  drop 
E  in  (2).  We  obtain 

( 3 )      ^r  "  '^r  ^  ^r  '  ^r  "  "  *^r^ 

where  y,,  ~  %J9^  -   ^r/^r*  "^^^  ^r  ®^®  *^®  gyro  -f re  quenc  ie s « 
Each  fluid  rotates  at  its  gyro  frequency,  CJ  ,      More  pre- 
cisely, each  fluid  translates  at  a  \iniform  velocity  u 
which  changes  its  direction  in  time.  The  velocity  com- 
ponents parallel  to  B  are  constant.  This  fluid  motion  is, 
of  course,  closely  related  to  the  motion  of  an  individual 
particle  in  a  uniform  magnetic  field,  namely  as  given  by 

A  non-zero  fluid  velocity,  u,  results  from  a  gas  of 

rotating  particles  when  the  phases  of  the  individual 

particles  are  not  completely  random. 

Now,  still  with  a   =  0,  we  set  B  =  0  but  keep  E. 
'  rs    ' 

Elimination  of  E  between  (1)  and  (2)  yields 

Pr^r  -  -  iT  ^  ^s  ^3   • 
Multiplying  by  \     and  siunmlng  yields 
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A  is   the  plasma  frequency.     Writing  J  =  Q.e  ,   we  have 

(E   =  -   j/iK.A 
Z  S.  u^  =  c 

The  basic  physical  effect  here  is  the  inertia  of  the 
electric  current,  i.e.,  of  the  particles  carrying  it. 
The  net  fluid  velocity  is  constant. 

As  In  the  case  of  gyro  oscillations,  the  motion  of  each 
fluid  is  a  translation  in  a  direction  which  varies  in  time. 
Considering  the  oscillating  part  only,  all  fluids  are  in 
phase.  The  polarization  is  arbitrary  as  compared  to  the 
previous  case  of  gyro  oscillations  In  which  the  polari- 
zation was  circular. 

If  the  oscillation  is  in  a  fixed  direction,  it  can 
be  considered  to  be  the  solution  to  the  problem  of  an 
oscillating  plane  condenser  discharge.   In  this  case,  the 
physical  interpretation  is  somewhat  more  intuitive j  there 
is  a  restoring  force,  E,  set  up  by  a  separation  of  charges 
from  the  fluid  into  the  condenser  plates. 
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The  conventional  formula  for  plasma  frequency  is 

where  m«  is  the  electron  mass.  This  is  an  approximation 
to  (i;)  obtained  by  considering  electron  mass  small  com- 
pared to  any  other  mass.   It  is  interesting  to  note  that 
there  is  no  separate  plasma  frequency  for  each  fluid  (as 
there  is  a  gyro  frequency);  there  is  a  single  plasma  fre- 
quency whose  value  is  dominated  by  the  particle  of 
smallest  mass. 

Next  we  turn  to  the  general  case  but  still  keep 
a   =  0.   It  is  easy  to  see  that  the  components  of  u 
and  E  parallel  to  B  and  perpendicular  to  B  separate  into 
two  uncoupled  sets  of  differential  equations.  The 
equations  for  the  parallel  components  are  identical  to 
the  plasma  oscillations  with  B  =  0  which  were  previously 
considered.   Therefore,  we  need  only  consider  the  perpen- 
dicular components  which  we  write 

t     " 

^r  ^  ^^r   ♦  ^r  ^* 

Eliminating  E  between  (1)  and  (2), 
(7)         ^r  =  -Tfl^s  ^s  ■"  Yr  ^r  ^  ^  • 


This  system  of  equations  can  be  written 

X     -  Y\         /  u' 

X 

-  7  - 


where 


Y  £  Y„  BD  5^^ 

D  =  ^ 


The  transformation 


y-i   x-i 


-x-i   y-i 


A  =  Y"-"-  X  +  X"-"- 


reduces  the  problem  to 


A  u'  =  0 


or 


Y"-^XAu'  =  (Y"-*-  X  +  1  I)  (Y"-'-  X  -  il)  u'  =  0,   1  =  5 


rs 


Prom 


Y"-*-  X  =  ~=r    (D^  6    +  Y  q  A) 
Y^BD      rs    'r  ^s'  '^ 

we  obtain  for  the  frequencies,  setting  D  =  iw  , 

(8)      \u){cj+0^)    5^g  -  Y^  q^  /k  I  =  0,  CJ^  =  -Yj.B 

Only  the  -  sign  need  be  taken  (changing  the  sign  of  cj 

merely  changes  the  sign  of  U)),  and  a  simple  computation 

yields  (dropping  zero  roots) 

1   n_  r^q 


K^  u)-a.^ 


or,  in  polynomial  form. 
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(9)    K   n  {  cj-  co^)   =^  YZ  r^  %  TTr  (  w  -  cjg) 

r=l  "^  r=l       s=-L 

n 
where  1  f  omits  s  =  r.  The  right  hand  side  is  a  poly- 

s=l 

nomial  since  the  sum  is  divisible  by  cj.   For  a  approach- 
ing zero  (low  density)  we  obtain  the  gyro  frequencies, 

(J  =  CJ  •   For  B  =  0  (  6J   =  0 )  we  obtain  the  plasma 
r  r  ^ 

frequency  and  the  root  CJ  =  0   with  multiplicity  2n  -  2, 
It  is  easy  to  expand  In  the  neighborhood  of  q  =  0  or 
(*>      =0.   In  the  first  case,  we  get 

(r)         "^r  %  % 

W    =  ^r  ^-^-='  •••  =  ^r  -KB  ^  ••• 

r 

and,  in  particular 

In  the  second  case,  the  double  root  •O'  splits, 

and  the  remaining  roots  are  first  order  in  Cu   and  given  by 

H.  Yj,  qp 

^  ^  =  0   . 


r"^  CO-(jp 

If  one  fluid  consists  of  electrons  and  all  others  are 

heavy  ions,  Y-i  ^^  Yp  •••  Y„  »  ^i  ^^  ^p  '••^n  »  "then 
there  are  two  modes  given  by 

u  {CJ  -cj-^)  =  n^ 

while  the  remaining  ones  are  of  the  order  of  magnitude 


of  U^    •..  CJ„, 

2       n 
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VJe  remark  that  the  general  solution  consists  of  a 
superposition  of  plasma  oscillations  parallel  to  B  and 
gyro-plasma  oscillations  perpendicular  to  B« 

It  is  important  to  realize  that,  in  the  general  case, 
the  quantities  CJ   and  /\     do  not  represent  natural 
frequencies  of  the  meditim,  but  parameters  from  which  the 
natural  frequencies  can  be  computed. 

To  treat  the  frictional  terms  in  any  sort  of  explicit 
form  it  is  necessary  to  restrict  the  discussion  to  the 
case  of  two  fluids.  We  have 

Pl^l  ~  "^1  (E  +  u^  X  B)  -  a  (u^  -  u^) 
^2^2  ~  ^2  (E  +  U2  X  B)  -  a  (U2  -  u^) 


(10) 


q^  +  q2  =  0 


Using 


'1  ^1 


P2  u^  =  p  u 


q^  u^  +  q^  U2  =  J  , 
we  introduce  u  and  J  instead  of  u^  and  Up, 


(11) 


ICE  =  -  J 


pu  =  J  X  B 


J  =  KA^  (E  +  u  X  B)  +  ^1^2  -^   ^2Pl  j^  _  _G^   J  ^ 

P  P1P2 


E  and  u  can  be  eliminated, 

(12)    y  =.a^j.}Ld.   ( j^)^  ^   ^lP2-*-r2Pl  ^^B  ,  _ae.  J  . 


P 
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Parallel  to  B  we  have 

(13)  J  -H(J^  J  +  ^2  J  =  0,    <<>d  =^ 

and  perpendicular  to  B,  introducing  Q.  =  J  x  B, 
(111)  y  +   CJ^^+    iXl^    ^   tJiU>2)Q.  +     ^  ^p    ^^  B  =  0 

As  before,  solutions  of  (13)  can  be  obtained  as  special 
cases  of  (li|),  taking  B  =  0  {  Cj^   =  CJ   =  O). 

The  frequencies  of  equation  (1^)  are  given  by 

(15)  O^  +  CJ  [iCO^  +  (U^  H-CO^)]  +^1^2  "^^  ""  °   » 

from  which  their  behavior  can  easily  be  seen.   In  the 
direction  parallel  to  B  we  have 

(16)  0^^  +  icj^ui  -  n^  =  0  . 

A  non-oscillatory  solution  of  (11)  has  J  =  0  and 
E  +  u  X  B  =  0.   The  general  solution  consists  of  an 
oscillatory  decay  towards  such  a  steady  solution  with 
damped  plasma  oscillations  parallel  to  B  and  damped  gyro- 
plasma  oscillations  perpendicular  to  B. 
2,      The  Debye  Length 

We  nov;  consider  the  static  fluid  distribution  which 
results  from  inserting  a  fixed  point  charge,  say  positive. 
The  tendency  will  be  to  repel  a  positive  fluid  and  attract 
a  negative  fluid.   For  simplicity,  we  linearize  and  ass\ame 
that  deviations  from  uniformity  are  small.  The  equilibrium 
equations  are 

-  11  - 


(17)  1        ''^"Y 

dlv  E   =  q/K 

and,   for  a  perfect   gas, 

(18)  Pp   =  v^kT^   =  q^kT/£^ 
so    that,    assxuning  T      is   constajit, 

r 
Smnming,   taking  the   divergence,    and  linearizing,   we   find 

Aq   =  q/d^ 
(20)  \  ^    _  ^2 


f         V 


1  =T  fli  =  V  

T^       ^  K-kT         ^    KkT  » 

a  r  r 

d   is   the  Debye   length.     The    solution  for   q  is 

-r/d 
(21)  "l  =  ='Vd-     • 

The  intuitive  interpretation  of  d  is  a  shielding  distance} 
beyond  a  distance  d  from  a  given  charge,  the  charge  is  not 
visible  (E  drops  to  zero  e:^ponentlally ).   The  fact  that  q 
becomes  infinite  as  r  >  0  violates  the  assumed  lineari- 
zation, but  essentially  the  same  results  can  be  seen  from 
the  full  nonlinear  equations. 

If  the  S,     are  equal  and  also  the  T   , 
r  r  ' 

2  ^  1^1^    RT 


d      2  '^   2 
(22)  )  ^^     ^ 

R  =  k/m 

where  m  is  the  electron  mass  and  the  approximation  (6) 
has  been  used.   The  connection  betx-ieen  plasma  frequency, 
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2 

Debye  length,  and  electron  speed  (RT'"*'v  )  would  seem  to 

be  intuitively  connected  with  the  charge  separation  inter- 
pretation of  plasma  oscillations. 

Another  property  of  d  comes  from 

I  ^^  ~    K       Kd 

(23)      \  3 

The  right  hand  side  can  be  interpreted  as  the  electro- 
static energy  between  a  charge  £  and  a  charge  N£  (the 
total  charge  in  a  Debye  cube  )  separated  by  a  distance  d. 
However,  the  potential  energy  assigned  to  a  particle 
because  of  interactions  with  its  neighbors,  i.e.,  as 
given  by  the  solution  of  (20),  has  the  order  of  magnitude 
e,  /Kd  since  particles  farther  than  d  are  ineffective,  and 
the  total  separated  charge  is  £•  We  conclude  that,  to 
have  a  perfect  gas  (g,  /icd  «  kT  ) ,  we  must  have  N  »  1; 
i.e.,  each  particle  must  be  in  the  potential  field  of 
many  others  I   At  normal  temperatiore  and  pressure  (which 
is  artificial  since  there  is  no  ionization)  N  -^2  x  lO"^, 

0 

while  at  a  pressure  of  1  mm  and  a  temperature  of  10,000  G, 
N  -^-lO,  which  is  a  fairly  good  perfect  gas. 

A  niomber  of  general  estimates  of  orders  of  magnitude 
can  be  made  using  the  lengths  L^-,  and  L^  defined  in  the 
last  chapter.  We  compute 

l2 


(21;)         a^  =  d^  2  ^  ^^  ^k 
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Frequently  £EL  has  the  same  order  of  magnitude  as  kT  if 

L  is  the  apparatus  size.   Furthermore,  d  is  usually  extremely 

-8 
small  (2  X  10   cm,  at  normal  temperature  and  pressure, 

3  X  lO"  cm,  at  10,000  C  and  1  mm.  pressure).  We  conclude 
that  q/vc  is  small,  i.e.,  the  medium  is  neutral  in  the 
sense  that  the  difference  between  the  number  of  positive 
and  negative  particles  is  small  compared  to  either.   Con- 
versely, if  this  is  not  so,  the  resulting  E  will  be 
enormous  compared  to  the  particle  energy.   Alternatively, 
a  small  domain  of  dimension  d  over  which  E  changes  appreci- 
ably (e.g.,  in  a  boundary  layer)  may  not  be  neutral. 

An  analogous  computation  can  be  made  for  the  magnetic 
field, 

H   ^  HB  ^  ^  _5. 

"  ^^^%   mvv2  %   % 

(25)   \        ^=fi=J 

mvv 

\   is  the  Larmor  radius,  namely,  the  radius  of  a  particle 
orbit  in  a  uniform  magnetic  field.   First  of  all,  since 
j/vev  cannot  be  larger  than  one  in  order  of  magnitude  if 
v  is  a  representative  electron  velocity,  we  find  that  Ltt 
cannot  be  smaller  than  5  in  order  of  magnitude.   If 
HB/mvv  "-^  HB/p  is  of  order  one,  the  magnetic  field 
cannot  alter  appreciably  over  a  distance  shorter  than 
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the  Larmor  radius.   Conversely,  If  Lj^  is  large  compared 
to  5,  the  relative  velocity  between  ions  and  electrons 
is  small  compared  to  thermal  velocities.   Comparing  with 
the  previous  example,  we  might  call  6  a  magnetic  shield- 
ing distance.   The  relation 

(26)  ^^JYr^ 

is  illuminating.   In  the  case  of  electrostatic  forces, 
d  separates  nearby,  i.e.,  collisional  type  forces  from 
long  range  forces.   In  the  case  of  magnetic  forces,  C 
accomplishes  the  same  purpose.   Talcing  HB/p  to  be  of 
order  one,  5 — ^—  d,  which,  combined  with  the  rough 
estimate  of  interparticle  magnetic  forces  as  being  on 
the  order  of  v/c  smaller  than  electrostatic  forces, 
implies  that  the  "collisional"  electrostatic  and  magnetic 
forces  have  equal  orders  of  magnitude. 
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